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Abstract 

With twisted boundary conditions on the quark fields, we study nucleon matrix elements of 
the axial current utilizing twisted heavy baryon chiral perturbation theory. One can explore the 
momentum transfer dependence of the axial form factors more easily than by using ordinary lattice 
quantized momenta alone. As examples, we derive expressions for the nucleon axial radius and 
pseudoscalar form factor. 
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I. INTRODUCTION 



Lattice QCD calculations will provide first principles determination of low-energy 
hadronic properties. This low-energy region of QCD is the regime of strong interactions, 
where perturbation theory fails and the quark and gluon degrees of freedom are confined 
into color neutral hadronic states. Advances in lattice QCD have largely come from ex- 
ploiting the features of numerically approximating the theory on a discrete lattice in a finite 
volume. As such, there is freedom to manipulate lattice quantities, so that the QCD answer 
is recovered in the continuum and infinite volume limits. 

A mutable feature of finite volume simulations is the choice of boundary conditions sat- 
isfied by the fields. Periodic boundary conditions are usually assumed as a matter of con- 
venience, as they lead to fields and hence observables which are single valued. Observables 
are determined, however, from S-matrix elements and hence only the action need be single 
valued for physical quantities to be well defined. The generator of any symmetry of the 
action can also be used to specify boundary conditions. In particular symmetries involving 
flavor generators lead to so-called twisted boundary conditions on the quark fields. The ap- 
pearance of twisted boundary conditions is not new. They have been considered in various 
contexts over the years [3, IE IE 0, IE IS IE IS E3] • Recently there has been renewed interest 
in utilizing twisted boundary conditions in lattice QCD simulations Q [13, EE El Q ■ This 
interest stems from the realization that the restriction to lattice quantized momenta can be 
bypassed by employing twisted quark fields. 

In this work, we consider the nucleon axial transition form factors. On a given lattice, the 
available momentum transfer is quantized, and on current lattices one has integer multiples of 
about a few hundred MeV. While chiral perturbation theory (xPT) predicts the momentum 
transfer dependence of these form factors, both the quark masses and lattice momenta must 
be brought down so that one enters the effective theory's range of applicability. Another 
difficulty encountered with a coarse sampling of the momentum transfer on a fixed lattice is 
the determination of radii and moments which are only accessible in the near-forward limit. 1 
These issues can be circumvented with twisted quarks. 

Before specializing to the case of the nucleon axial matrix elements below, we begin 
with a more general observation about single particle matrix elements of flavor-changing 
operators. 2 These matrix elements have the form (H'(P')\0\H(P)). Here the operator O 
changes the flavor composition of the hadron H into that of H', and we assume at least 
one of the states is not flavor neutral. The above matrix element can be decomposed into 
various form factors which depend on the momentum transfer. Let the size of the lattice be 
L in each of the spatial directions. On the lattice, the momentum transfer q = P' — P is 
quantized q = 27cn/L, for n G Z 3 , when one imposes periodic boundary conditions on the 



1 To sample smaller values of the momentum transfer, one could work with lattices that are much longer 
in only one spatial direction. With such lattices, however, working with sufficiently small lattice spacings 
and quark masses is well beyond the reach of current simulations. 

2 A similar observation can also be made for particle to vacuum transition matrix elements. For a generic 
pseudoscalar H, we have schematically (0\A\H(P)) = ifnP, where A. is the spatial part of the axial- 
vector current. With twisted boundary conditions and for a flavor non-diagonal state H , one can access 
the decay constant fn at zero lattice momentum. While one could have easily deduced the decay constant 
by working in the rest frame and using the time-component of the current, the above observation provides 
a novel check on using twisted quarks. 
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quark fields. To smoothly access the momentum transfer dependence of the transition form 
factors, one is forced to vary the lattice size and hence generate new gauge configurations. 3 
Imposing twisted boundary conditions on the valence quark fields, such that different flavors 
are independently rotated in phase at the boundary, one finds an induced momentum transfer 
q = 2nn/L + 50/ L, where 50 is the difference in twist angles of the flavors changed. Thus 
with twisted boundary conditions on the valence quarks, the momentum transfer can be 
continuously varied and one can access off diagonal matrix elements, even at zero lattice 
momentum transfer n = 0. For such matrix elements, one can work in the near forward 
limit without utilizing a box with one very long side, and without performing a momentum 
extrapolation. 

In the following section (Sec. HI]), we focus on the nucleon axial matrix elements. Here we 
extend our momentum transfer observation above to the case of the neutron-proton axial 
correlation function. Next in Sec. 11111 we develop partially twisted baryon for lattice 
calculations in which the valence (and ghost) quarks satisfy twisted boundary conditions, 
while the sea quarks remain periodic. This effective theory can be used at finite volume 
to ascertain the corrections to nucleon matrix elements due to twisting. As an application, 
we calculate the axial radius and pseudoscalar form factor of the nucleon in Sec. IIV[ and 
summarize our findings in Sec. 

II. NUCLEON AXIAL TRANSITION 

Let us now focus specifically on the nucleon axial matrix elements and make our introduc- 
tory observation more concrete. For simplicity, here and below we work in the isospin limit, 
and consider the neutron and proton to be degenerate. We also frame our discussion using 
continuum operators. The relevant flavor-changing operator for the nucleon axial transition 
is O ll (x) = u{x) r ) ^^d{x) . The Euclidean time, axial correlation function has the form 

(V (t)O»(t'Wrm = Y. e ~ lPX ' (0\V(x,t)O,(x',t'W(0,0)\0), (1) 

x,x' 

where the sums are over all lattice sites labeled by x, while M{x, t) and V(x, t) are interpo- 
lating fields for the neutron and proton, respectively. On the left-hand side of the equation, 
the subscripts on the interpolating fields denote the spatial momentum. In writing this 
correlation function, we have assumed that t > t' > 0, and for simplicity have chosen the 
final state to be at rest. As is standard, the ground state nucleons will be filtered out in the 
limit of large Euclidean time separation, t ^> t' ^> 0, and the nucleon axial matrix element 
can be isolated by taking a ratio of the three-point and two-point functions weighted with 
the appropriate kinematic factor. 

Due to the periodicity of the fields, the initial-state momentum is quantized P = 2-n-n/L, 
and the available momentum transfer probed on a given set of gauge configurations is hence 
coarsely grained for current lattice sizes. To circumvent this restriction, we replace the 



3 In principle, one could use the dynamically generated quark mass dependence of the initial and final state 
hadron masses to explore the momentum transfer dependence of such form factors, assuming that these 
hadrons are non-degenerate. This does not give one direct control over the momentum transfer probed, 
moreover, the procedure is complicated by the additional need for chiral extrapolation. 
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interpolating fields built from periodic quark fields with new interpolating fields formed 
from twisted quark fields. We shall use the same notation for these fields used above: 
Af(x, t) and V(x, t). The precise form of the twisting will be spelled out later in Sec. IIHI It 
suffices to note that each valence quark flavor independently suffers a discontinuous change 
of phase at the lattice boundary, and we denote 6 U and 6 d as the twisting angles for u and 
d quarks, respectively. 

To implement the twisted boundary conditions, one uses modified fields that are periodic 
but coupled to a £7(1) gauge field. For the quarks, one calculates their propagators in this 
background gauge field |l2j |. where each flavor's twisting angle acts as an induced charge. 
For the baryons, we can formally separate off the kinematic effects of twisting by defining 
the periodic interpolating fields 

M{x,t) = e- iB «- x N{x,t), 

V(x,t) = e- iBvX V(x,t), (2) 

where B N = (0 U + 20 d )/L and B v = (20" + d )/L are the effective £7(1) charges of the 
neutron and proton, respectively. Using these periodic fields, one then calculates the various 
Wick contractions in terms of modified quark propagators. 

In order to determine the axial correlation function, one must perform the Wick con- 
tractions when the axial operator is inserted 4 . To do so, one inserts the similarly modified 
axial- vector operator O^x) which is defined in terms of periodic quark fields with effective 
£7(1) charges. This operator is related to the twisted quark axial- vector operator via 

6,{x) = e- i{ * ed - eu > x l L 0,{x). (3) 

The calculation of Wick contractions then proceeds without any extraneous position- 
dependent phase factors. To obtain the correlation function using twisted fields on the 
lattice, one hence determines 

e~ lPx ' (0\V(x, t)6,(x', t')Jf(0, 0)|0) 

x,x' 

= e - iBv < x - x,) e- l{p+B ^- x ' (0\V(x, t)0„(x', t'W(0, 0)|0) 

x,x' 

= (v Bv {t)o tl {t')jr P+BM {Q)), (4) 

where in the subsequent steps we have rewritten the lattice correlation function to expose 
that the initial and final states have been boosted. The momentum transfered to the fi- 
nal state is effectively q = (9 U — 6 d — 2irn)/L and remains non-zero even at zero lattice 
momentum, i.e. when n = 0. 

While the momentum transfer in such matrix elements can varied continuously, twisting 
produces long-range flavor symmetry breaking interactions that modify the physics we seek 
to explore. On the lattice, twisting thus introduces modified finite volume corrections that 
can be determined using chiral effective theories 0|. To this end, we develop partially 
twisted heavy baryon chiral perturbation theory and apply it to the nucleon axial matrix 
elements. We ultimately address two simple examples that do not have sizable finite volume 
corrections. 



4 Notice there are no self-contractions of the operator (x) . 
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III. PARTIALLY TWISTED BARYON CHIRAL PERTURBATION THEORY 



To address the consequences of twisting in lattice calculations of baryon properties, we 
construct the underlying effective theory in the baryon sector. First we detail the partially 
twisted boundary conditions employed and then proceed to include these effects in heavy 
baryon xPT. The quark part of the partially quenched QCD Lagrangian is 

6 

£= J2Q j ( i] P- m Qh k Qk- (5) 

j,k=l 

The six quark fields transform in the fundamental representation of the graded S77(4|2) 
group and appear in the vector Q T = (u,d,j,l,u,d). In addition to the u and d 
quarks, we have added ghost quarks u and d, which cancel the closed valence loops, and 
two sea quarks j and I. In the isospin limit, the quark mass matrix of SU(A\2) reads 
tuq = dmg(m u ,m u ,mj,mj,rn u ,m u ), so that QCD is recovered in the limit rrij — > m u . We 
require that the quark fields satisfy twisted boundary conditions, namely 

Q(x + Le r ) = exp «T a ) Q(x), (6) 

where e r is a unit vector in the r th spatial direction and the block diagonal form of the 
supermatrices T a is 

T a = diag(T a ,0,T a ). (7) 

Here T a are the elements of the U(2) algebra. In the isospin limit, any generator can be 
chosen for the twists, although we choose to preserve electric charge conservation and ac- 
cordingly restrict T a to the Cartan subalgebra. Notice in Eq. (J7J) the sea quarks remain 
periodic at the boundary. Consequently the twist angles can be changed without necessitat- 
ing the generation of new gauge configurations and the fermionic determinant, which arises 
solely from the sea sector, is not affected by the twisting. 

Redefining the quark fields as Q(x) = V'{x)Q(x), where V(x) = exp(iO a ■ xT a /L), we 
can write the partially quenched QCD Lagrangian as 

c= Y,Q 3 (^- m o) (8) 

j,k=i x 3 

where all Q fields satisfy^ periodic boundary conditions, and the effect of twisting has the 
form of a gauge field: = + iB^, where = (0, 6 a T a /L). It will be easier to 
treat the twisting in the flavor basis of the valence and ghost sectors rather than in the 
generator basis, thus we write 6 a T a = diag(6 u ,6 d ,0,0 J 6 u J 6 d ), and similarly for B^ = 
diag(^,^,0,0,S«,^). 

In the meson sector of partially quenched xPT (PQxPT) [3, [l?], EH EM El; the coset 
field E, which satisfies twisted boundary conditions, can be traded in for the field E defined 
by E(x) = (x)T,(x)V (x) , which is periodic at the boundary jl4]. In terms of this field, 
the Lagrangian of PQxPT appears as 

£ = ^-str (D^ED^ + A str (m ] Q t + S+mg) . (9) 

The action of the covariant derivative _D M is specified by Z) M E = <9^E + i[B^, E]. 
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To include baryons into PQxPT, one uses rank three flavor tensors [2l(. The spin-l 
baryons are described by the 70- dimensional supermultiplet B l ^ k , while the spin-| baryons 
are described by the 44- dimensional supermultiplet T^ k |22|. The baryon flavor tensors are, 
however, twisted at the boundary of the lattice. In the r th spatial direction, both tensors 
satisfy boundary conditions of the form 



B ijk (x + e r L) = (-yAvj+v^+iVi'+Vi'Xm+w) tffk 

\ / ii' \ J jj' V / kk' 



X) 



e< r e^M U^ Ta ) B ijk (x), (10) 
/ a v / jj v J kk 

where the 77's are the grading factors of SU(A\2). In the last line, we used the diagonality of 
the Cartan generators. Thus we define new tensors B^ k and both having the form 

B ijk (x) = Vl{x)V^{x)Vl k {x)B ljk {x). (11) 

These baryon fields satisfy periodic boundary conditions and their free Lagrangian has the 
form 



£ = i (Bv ■ VB) + 2a ( M Q) (bBM+^J + 2f3fj Q) (t3M + B) + 2a { ^ Q) (Bb) str (M 

-i (7> ■ Vf^ + A (t^) + 2 7 ^ Q) (f,M + f^ - 2a™ (f^ str (M+ 

(12) 

where the mass operator is defined by M.+ = | (i^rriQ^ + irriQ^j, with £ = and the 
covariant derivative acts on B and 7^ fields in the same manner 

[V,J3{x)Y ik = d^& jk (x) + i{Bl + B^ + B k )B ljk (x). (13) 

The connection of the low-energy constants appearing in Eq. ()12j) to those of SU (2) %PT is 
described in [22]. 

With twisted boundary conditions, the leading order PQxPT interaction Lagrangian 
between the baryons and mesons now appears as 

C = 2a (fiS^BA^ + 2(3 (Is^A^ + 2H {f v S^A^f v ^ + (f v A v Bj + (&A V % 

(14) 

where the twisted axial- vector meson field is defined by A^ = | (^D^^ — ^D^^j. The 

familiar low-energy constants of SU(2) can be identified as follows [22|: g& — |a — 
gAN = —C, and #aa = Ti- Notice there is an extra free parameter in the PQxPT interaction 
Lagrangian compared to that of 



IV. NUCLEON AXIAL RADIUS AND PSEUDOSCALAR FORM FACTOR 

The axial current matrix elements of baryons have been studied extensively in chiral 
effective theories and many investigations have now been tailored in an attempt to describe 



the numerical approximations in lattice QCD calculations, see |22j, |23j, |24|, |25|, |26|, |27|, |2 
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29L 30l 3l| . Determination of the nucleon axial char 



37, 38, 39, 40, 41, W2l W 



44 l45j , but have not been investigated on 



investigations, for example |32j, |33|, |3J, |35|, |3J 
axial form factors have been calculated in ;\PT 
the lattice. Here we consider the nucleon axial form factors in partially twisted xPT. We 
detail how the axial radius and pseudoscalar form factor can be accessed, even at zero lattice 
momentum. 



ge has been t he g oal of numerous lattice 

The nucleon 



In partially quenched QCD, the iso- vector axial- vector current is defined by J 5 



Q7 M 7 5 (r a /2)Q. The choice of supermatrices ~T a is not unique [46j, even when one imposes 
the condition strr a = 0. One should choose a form of the supermatrices that maintains 
the cancellation of valence and ghost quark loops with an operator insertion ji?}- For the 
flavor changing contributions we consider below, however, these operator self-contractions 
automatically vanish. Thus regardless of the form chosen in extending t°, the only discon- 
nected quark contributions are those from the gauge configurations. For our calculation we 
require the action of in only the valence sector, and specify the upper 2x2 block of 
T a to be the usual Pauli isospin matrices r a = (1,t). Henceforth we restrict our attention 
to the operator = j}^ + i = u / y fJ- 'j 5 d. The no n- vanishing nucleon matrix element of 
this operator is the neutron to proton axial transition, which we decompose into two form 
factors 

q^q-S 



{p(P>)\J^\n(P))=u(P>) 



2S»G A (q 2 ) + 



(2M 2 



N , 



;G P (q 2 



u(P), 



(15) 



where q^ = (P' — P) M is the four-momentum transfer. Above G A (q 2 ) is the axial form factor 
and Gp(q 2 ) is the induced pseudoscalar form factor. The nucleon axial charge is G A (0), 
and in the chiral limit G A (0) = g A • Kinematically one can probe zero momentum transfer 
because in the isospin limit the proton and neutron are degenerate on the lattice, as the 
electromagnetic mass difference is absent. We define the axial radius < r\ > from the small 
momentum transfer expansion of the axial form factor 



G A {q 2 ) = G A ^) + q -<r\ >+... 
o 



(16) 



Clearly to access the pseudoscalar form factor and axial radius, one needs momentum trans- 
fer between the initial and final states. 



In partially twisted xPT, the form of the leading-order baryon axial current J 5 is 
J 5 +M = 2a (BS^BT + ^j + 2f3 (bS^t+B) + 2H (t u S^t + T^ 
+ a/|c Kt^t+b) + (BT + f^ 



(17) 



where the full chiral structure of the axial- vector current can be obtained under the replace- 
ment r + — > \ (<f^ + f + f r + f * J . To find contributions to the axial radius and pseudoscalar 

form factor, we must consider tree-level contributions of which there is the pion pole term 
and local interaction terms. The latter are derived conveniently from the Lagrangian. Let 
be the field-strength tensor of the external axial field. In untwisted xPT, one has the 
dimension-six interaction term contained in the Lagrangian 



£ = ^ (NS,r + N) d u F : 



A 2 

x 



hi/ 

A ' 



(18) 
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that leads to a contribution 5Jg to the axial current at next-to- leading order 



5,n 



2ua 



[d»d u (NS„t + N) - d 2 (NS^r+N)] . 



In partially twisted ;^PT the analogous next-to-leading order axial current reads 



(19) 



6J+» = ±l2n, 



b^b u (bs v t + b) - b 2 (bs^t+b 




(20) 



The relation of tia to the PQ^PT parameters can be found from matching, tia = \n a — \tip. 

Also, the terms of the axial current in Eq. (j!7|) can be combined with those of the 
interaction Lagrangian in Eq. (fTf)) to generate loop contributions to the axial radius and 
pseudoscalar form factor. These two observables, however, are rather special in heavy baryon 
%PT, because tree-level terms dominate over pion loop contributions. This is because the 
axial current is inserted on the baryon lines and generates terms ~ q 2 /M"^, rather than terms 
~ q 2 /m 2 b present in electromagnetic form factors, for example. Thus these quantities are 
largely insensitive to the long-range effects introduced by twisting. Subsequently twisting 
can be used to produce momentum transfer between nucleon states without sizable finite 
volume corrections. In calculating these quantities, we find 



<r\> 



1 + 



G P {q 2 ) = (2M 



Nj 



9a 



q 2 — ml 



2ua 

A 2 

x 



n , ^| q*_ ^ 



N 



(21) 
(22) 



where the recoil corrections (oc M N ) can be determined at the one-loop level [2a], while 



the equally sizable chiral corrections (oc A~ 2 ) enter at two loops [48]. Above is an 
abbreviation for the meson masses of PQ^PT. Because these PQ^PT results are determined 
at tree-level, sea quark masses do not appear and the dependence on the parameters n a and 



np enters only through the combination 



-,n r 



\n 



ha- In this way, the partially quenched 



results maintain the Adler-Dothan relation for the pseudoscalar form factor [49j . 

In writing Eqs. ()15j) . (|2T| . and (|22|) . we have used a compact notation for the momenta 
P and P' which have the form 



(2ixn 



'M 2 N +[--^ + B N 



( 2iim 



B 



v 



2nn „ 
, — + B M 



27rm „ 



(23) 



(24) 



as well as for the momentum transfer, which is g M ~ (0, P' — P) in the non-relativistic 



S 



limit. 5 We have chosen the above notation to be consistent with Sec. |H] Finally let us spell 
out our results and the resulting form of Eq. ()15|) at zero lattice momentum, i.e. m = n = 0. 
To this end we define B^ + = (B-p — BjsfY = (B u — B d )^, which is the induced U(l) charge 
of the 7r + field due to twisting. At zero lattice momentum, we have 



<p(0)|J 5 + >(0)> = u(0U2S» 



G A {0) 



1 



-B^ + B w + ■ S 



9a 1 2 

Bl + + ml + 3 <rA > 



>u(0) 



(25) 



where we have used < r\ > as a replacement for the quantity in Eq. (|2*T|) . This form makes 
transparent how varying the twisting parameters of the quarks can lead to a determination 
of the axial radius and pseudoscalar form factor. We remark that these determinations 
suffer additional finite volume corrections but these do not occur until higher orders in the 
heavy baryon and chiral expansions. The simplicity of these results, moreover, can be used 
as a consistency check on the chiral extrapolation. For example, if one sees valence or sea 
quark mass dependence of the axial radius greater than ~ m^/Mfj x 100%, then one is 
certainly not in the chiral regime. The same is true of the sea quark mass dependence of 
the pseudoscalar form factor, which has, however, dramatic behavior as a function of the 
valence quark masses. 



V. SUMMARY 

For hadronic matrix elements of flavor changing operators, twisted valence quarks produce 
momentum transfer between the initial and final states. Thus such form factors' momentum 
transfer dependence can be explored without the additional generation of gauge configu- 
rations. Moreover the near forward limit can be explored to extract radii, etc., without 
necessitating a box with a very long side and subsequent momentum extrapolation. There 
are, however, modified finite volume corrections to these form factors arising from twisted 
boundary conditions. For large enough volumes and below particle production thresholds, 
these corrections are exponentially small. We develop partially twisted baryon ;\PT as a 
means to address the extrapolation of lattice QCD calculations employing twisted valence 
quarks. We pursue two simple examples, the nucleon axial radius and pseudoscalar form 
factor. These quantities do not receive finite volume corrections until higher orders in the 
chiral expansion, thus providing a clean test case for the utilization of the induced mo- 
mentum transfer from twisted quarks. Moreover, the applicability of in the baryon 
sector can be tested due to the simple quark mass dependence predicted at leading order for 
these observables. Lastly confrontation with experimental data (which is ample for these 
observables) presents an essential test of lattice techniques. 



5 Here we have neglected the B dependence of the proton and neutron masses that is dynamically generated 
from finite volume effects. These effects lead to M n — M p ^ even in the isospin limit, but the difference 
introduced is beyond the order we are working. 
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